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ABSTRACT 

We determine the critical point and the ratios l3/v and of critical ex- 
ponents of the deconfinement transition in SU{2) gauge theory by applying the 
X^-method to Monte Carlo data of the modulus and the square of the Polyakov 
loop. With the same technique we find from the Binder cumulant Qr its universal 
value at the critical point in the ther mo dynamical limit to —1.403(16) and for the 
next-to-leading exponent lj = 1 ± 0.1. From the derivatives of the Polyakov loop 
dependent quantities we estimate then The result from the derivative of Qr 
is l/u = 0.63 ± 0.01, in complete agreement with that of the 3d Ising model. 
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1. Introduction 

The investigation of the critical properties of any physical theory on a finite 
lattice is a non-trivial problem. This is so, because the singularity of the infinite 
volume theory at a phase transition is modified due to the finite volume into 
a rapid but nonsingular variation of the order parameter. The location of the 
critical point becomes then unclear and the infinite volume form of observables 
with critical behaviour like the order parameter, the susceptibility and higher 
cumulants cannot be used for the determination of the critical exponents. The 
observed finite size effects are on the other hand controlled by these exponents, 
as described by finite size scaling (FSS) theory. Simulations on lattices of varying 
size allow then to calculate the critical parameters of the theory. 

Recently, a new method, the x^-method, was proposed [1] and applied to the 
expectation value of the square of the magnetization. Though this quantity is not 
peaked at the transition, it was shown that one may deduce from its FSS behaviour 
the precise asymptotic critical point and the ratio of the critical exponents 7 and 
u in a very straightforward manner. 

In this letter we want to extend the method to the modulus and Binder's 
fourth-order cumulant Qr [2] of the magnetization. The intention is to find the 
ratio P/u of exponents and to check the consistency of the results for the critical 
point from the different variables. As it turns out it is possible to obtain still 
further information on universal quantities and by examining the derivatives of the 
variables at the critical point one may even determine the inverse of the exponent 
u of the correlation length with - at least for SU (2) - unreached precision. 

2. Finite size scaling theory and the x^-method 

Let us briefly review that part of FSS theory which is relevant for our con- 
siderations. In the neighbourhood of the critical temperature Tc of a second order 
phase transition one expects in the limit of large A^^-, the characteristic length 
scale of the system, for the magnetization or order parameter (denoted by (L) 
here) that 

(L)~(T-Te)^, for T^T+. (1) 

The behaviour near to Tg of the susceptibility x ^■nd the correlation length ^ in 
the large No- limit is expected to be 

X-IT-Tep^ (2) 
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e-|T-Te|-^ (3) 

A quantitative analysis of the finite lattice modifications to these formulae becomes 
possible by using renormalization group theory. In this framework it has been 
shown [3] that the singular part of the free energy density has the following form 

fs{x,h,N,) = N-''Qf^{gTNy'^,gf,N(^+^y\g,Ny^) . (4) 

The scaling function Qf^ depends on the temperature T and the external field 
strength h in terms of a thermal and a magnetic scaling field 

grp = ctx + 0{xh, x^) , (5) 

Qh = Chh + 0{xh, h^) , (6) 
which are independent of N^j and where x is the reduced temperature 

T — T 

c,oo 

Here the action contains a further symmetry breaking term h ■ N^L. Also 
additional irrelevant scaling fields gi with negative exponents may be present. 

The order parameter (L), the susceptibility x the cumulant gr are ob- 
tained from fs by taking derivatives with respect to h at h — 0. The general form 
of the scaling relations derived in this way is 

0{x,N,) = N^/^ ■ QoigrN^^giNy^) • (8) 

Here O is {L),x ^^^d gr with (p = — 7 and 0.) Taking into account only the 
largest irrelevant exponent yi = —uj and expanding the scaling function Qo to 
first order at a; = results in the following equation 

0{x, N,) = NP'^[co + (ci + C2N-^)xNl'^ + csiV"-] . (9) 

Standard FSS methods are based on the evaluation of Eq. (8) and/or Eq. (9) in 
the neighbourhood of the infinite volume critical coupling. 

The x^-method [1] utilizes Eq. (9) only at a; = 

0{x = 0, N,) = NPj^[co + csiV-'"] , (10) 



2 



and determines the critical point from the leading A^cr-behaviour of Eq. (10). This 
is motivated by the fact, that for x 7^ the A^^r -behaviour is drastically changed 
due to the presence of xNy^-terms ( for SU{2) e.g. 1/u ^ 1.6 [4]). The critical 
point is then defined as that point where a fit to the leading A'^^r-behaviour of 
Eq. (10) has the least minimal x^. 

We have two different forms of fits. If the exponent p 7^ the leading A^^^- 
behaviour is given only by the first term in Eq. (10). Taking then the logarithm 
we find 

lnO = lnco + ^InAT^ , (11) 

i.e. we have a linear dependence on InA^g. with slope p/v. Linear fits to In O as a 
function of In A^^- give then also the exponent ratio. A different form of fit is to be 
used in the case p = 0, where the leading A^o-behaviour is 

O = co + C3N-^. (12) 



It is interesting to consider also derivatives of the observables with respect to 
X. Since the scaling functions Qo are depending on xNa^^ we obtain 

^ix,N^) = Ni'+^y'^ . Q'oi^N'Jn , (13) 

where the prime denotes the derivative with respect to the argument; the de- 
pendence on the irrelevant fields is not explicitly shown, but present. Since the 
derivative of the scaling function is again a scaling function, derivatives of ob- 
servables O may be used to find the value of (1 -|- p)/iy from fits to the leading 
A^o--behaviour of 

^(x = 0, N,) = Ar(i+'')/-[ci + . (14) 

Derivatives may also be used to define new direct scaling functions by 

U = = UixN^n ■ (15) 

If for any finite value of A^^r the scaling function Qo and its derivative Q'q are finite 
in a neighbourhood of a; = - which is the normal case - we have immediately 

U{0) = 0. (16) 
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In the thermodynamic hmit, i.e. for A^^r oo, we get, however a different result. 
From Eqs. 1-3, or in general from 

Ooo = cox-P for , (17) 

we find 

U^{0) = -p. (16) 

An attempt to determine p at a; = from infinite volume formulae will therefore 
usually fail on finite lattices. Yet at small, but finite, positive x and large A^o-) i-^- 
for large arguments xN^^^ , U will approach —p. 



3. SU(2) gauge theory at finite temperature 

Let us consider SU (2) gauge theory on x N^- lattices using the standard 
Wilson action 

^ p 

where Up is the product of link operators around a plaquette. The number of 
lattice points in the space (time) direction A^cr(T) ^-nd the lattice spacing a fix the 
volume and temperature as 

V = {N^af, T = l/{Nra) . (19) 

On an infinite volume lattice the order parameter or magnetization for the 
deconfinement transition is the expectation value of the Polyakov loop 



^ AT, 

L(x) = -Tr n t/.,x;4 , (20) 



T = l 



or else, that of its lattice average 

^ = i3E^W' (21) 

where Ux-4 are the SU{2) link matrices at four-position x in time direction. 

Since, due to system flips between the two ordered states on finite lattices 
the expectation value (L) is always zero, the true susceptibility 

X ^ KiiL') - (Lf) , (22) 
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reduces there to 



(23) 



The quantity Xv is monotonicaUy rising as a function of the couphng A/ g"^ or the 
temperature T. Below the critical point 

X. = X , (24) 

and one expects therefore Xv to have the same FSS behaviour as the true suscep- 
tibility for T<T(.. In rcf. [1] it was shown that this is indeed the case. 

In the following we want to extend the analysis of ref. [1] to the remaining 
magnetization dependent observables and their derivatives. In particular, we want 
to take the expectation value of the modulus of the lattice average, as a kind 

of "finite lattice order parameter", hoping as usual [4,5] that its FSS behaviour is 
controlled by the exponent (3 of the infinite volume order parameter. In addition 
we want to investigate the behaviour of the normalized fourth cumulant of the 
magnetization 

/r4\ 

9r = -^—^ - 3 , (25) 

which is directly a scaling function and corresponds, up to a constant, to the 
renormalized coupling defined for infinite systems [6]. 

4. Numerical results and discussion 

The SU{2) Monte Carlo data, which wc want to use here were computed 
on A*"^ X Nj- lattices with A^^- = 8, 12, 18, 26 and A^^ = 4 and have already been 
reported on in refs. [5,7] and [1]. Since in the latter reference the critical point 
was determined to = 2.2988(1), we calculated in addition the still missing 

results at this point for A^^- = 8, 12 and 18 and A^^- = 4 with high statistics. They 
are presented in Table 1. Proceeding in the same way as in ref. [1], we evaluated 
then the data in the very close vicinity ( 2.298 < A/g'^ < 2.300) of the transition 
with the density of states method (DSM) [8]. The newly calculated MC data were 
in complete agreement with the former DSM interpolation and introduced only a 
tiny shift to the interpolation itself. As usual, the errors were determined with the 
Jackknife method. 
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(1^1) 




gr 


8 


4 


0.14343(18) 


0.025723(54) 


-1.3532(20) 


12 


4 


0.11591(22) 


0.016749(48) 


-1.3695(37) 


18 


4 


0.09370(88) 


0.01091(15) 


-1.386(16) 


26 


4 


0.07649(75) 


0.00733(10) 


-1.365(17) 



Table 1 

MC results at = 2.2988. 

At each 4/^^-value we have made hnear x^"fits to our DSM interpolation results. 
For (|L|) and Xv the form of Eq. (11) was used, i.e. ln(|L|) and Inx^, were fitted 
as a function of InAr^.. Correspondingly, Eq. (12) was utilized for fixed values of uj 
near to 1 to perform linear fits of ctS ct function of A^^ . In Fig. 1 we compare 
the obtained minimal x^/^f for the three quantities. Here the number of degrees 
of freedom, Nf, is 2, and uj was fixed to 1. 

We observe, that all three observables give consistent results in several re- 
spects : the positions of the minima of the x^^^/Nf-cmYes are nearly coinciding - 
their differences give us an error estimate on the critical coupling; the curves have 
the same shape and show a very sharp increase, symmetric around the respective 
minimum; the size of Xmin/-^/ close to the minimum is of order one, which indi- 
cates that our error estimates are reasonable. Finally, the minima themselves are 
very close to zero, which means that the finite size behaviours of the observables 
at the critical point are indeed described by the formulae, Eqs. 11 and 12. In Table 
2, we list the minimum positions found. 



o 


• 


UJ 


m) 


2.29895 




Xv 


2.29890 




9r 


2.29905 


0.9, 1.0 




2.29900 


1.1 



Table 2 

Minimum positions of Xmin/-^f- 
We conclude from their variations, that the critical point is at 

^ 2.29895(10) . (26) 
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Inside error bars, this value is compatible with the one of ref. [1], which was 
determined from the finite volume susceptibility Xv only. The shift in ^/ g'^iniXv) 
due to the increased statistics is 0.0001. As can be seen from Table 2 a variation 
of uj in the range 0.9-1.1 has nearly no effect on the position of the minimum. 
Increasing uj further shifts the minimum point to smaller 4/ (/^-values, the width 
of the x^j„/A^/-curve becomes smaller, but at the same time the minimim value 
increases slightly. Decreasing u) leads to opposite effects: shift of 4:/g'^-^ to higher 
values, widening of the x^j^/A^/-curve, decrease of the minimum. A value of u = 
1.0 or slightly higher seems reasonable, because then the width of the Xmin/-^f~ 
curve is the same as that for {\L\). Also the minimum position is then in good 
agreement with that of the other observables. In ref. [7] , where the critical points 
for Nr — 4 and 6 were determined from the intersection points of the (/^.-curves 
for different N^j, a value of c<; = 0.9 was found for A^,- — 6, whereas for Nj- — 4 the 
contribution of irrelevant scaling fields was below visibility. With our approach 
we confirm the existence of these fields also for N.,- = 4 with an exponent uj of 
comparable size. In addition, a; ~ 1.0 is the value found in the 3d Ising model 
[4], which belongs to the same universality class as 3+1 dimensional SU{2) gauge 
theory [9]. 

Each of our three fits yields one universal quantity: l3/v and ^/v are obtained 
from the (|L|) and Xv fits; the constant cq from the Qr fit is the universal value 
of gr at the critical point for N„ —>■ oo. Since all universal quantities should be 
equal for theories of the same universality class, we compare our results to those 
of the Sd Ising model [4]. This is shown in Fig. 2. Obviously we find excellent 
agreement between the two theories. In Table 3 we have listed the corresponding 
numbers at ^/ gl^^o- Their errors come from two sources: the uncertainty in 4/gf^ 
and the fit errors at fixed A/g"^. The value of g"^ was estimated in ref. [7] from 
Binder's [2] "cumulant crossing" method to -1.38(5). Our fit leads to a lower 
result -1.403(16), which is in closer agreement with the Ising value -1.41. This 
comes about, because we include a correction to scaling term in the fit (csA^"'^), 
where C3 turns out to be positive. 

As a final application of FSS methods we have investigated the derivatives of 
our observables with respect to the coupling 4/fif^. Since the reduced temperature 
X may be approximated by 

4/^2 _ 4/ 2 
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we have proportionality between 80/ dx and 90/9(4/(7^). As a consequence 
Eq. (14) may be used to find (1 + p)/^ from a fit of the type of Eq. (11) at 
the critical point. 



Source 




SU{2) 


Ising [4] 


(1^1) 




0.525(8) 


0.518(7) 


^(1^1) 


(1 - 


1.085(14) 


1.072(7) 






1.610(16) 


1.590(2) 




V 


0.621(6) 


0.6289(8) 






0.326(8) 


0.3258(44) 


Xv 




1.944(13) 


1.970(11) 


DXv 


(1 + 7)A 


3.555(15) 


3.560(11) 






1.611(20) 


1.590(2) 




V 


0.621(8) 


0.6289(8) 




7 


1.207(24) 


1.239(7) 






2.994(21) 


3.006(18) 


9r 


-9^ 


1.403(16) 


1.41 


Dgr 


Ijv 


1.587(27) 


1.590(2) 


(a; = l) 


V 


0.630(11) 


0.6289(8) 



Table 3 

Results from the x^-method at = 2.29895(10). 

Here, DO denotes the derivative dO / d{A/ g^). 

Whereas it is relatively easy to calculate the derivatives from spline interpolations 
to the DSM results, it seems more difficult to estimate their errors. To achieve 
this, we applied the Jackknifc method to those spline derivatives, which were 
obtained from the individual Jackknife block results. In Fig. 3 we show the slope 
parameters from the derivative fits in the vicinity of the critical point. Again we 
find a surprising agreement with the 3(i Ising model. In Table 3 the corresponding 
numbers are given as well as the resulting l/z/(z/) and exponent values. At the 
critical point the minimal /Nf was of the order 3 — 4 for the derivatives of {\L\) 
and Xv and about 1 for the derivative of g^. We think therefore, that the result 

V = 0.630(11) . (28) 

From the derivative of gr is the more reliable one, though all results are compatible 
with each other. We have checked also the hyperscaling relation 

-i/u + 2(3/u . (29) 
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As can be seen from Table 3 our results are perfectly consistent with Eq. (29). 

In summarizing the following comments are in order. We have shown, that 
also for the more complex (as compared to the Ising model) SU{2) gauge theory 
it is possible to determine the universal quantities related to the deconfinement 
transition with high precision from simulations on finite lattices. The special FSS 
method which we applied here, the x^-method, proved to be simple and accurate 
and showed consistency in the determination of the critical point from different 
observables. An indispensible part of the method is the DSM interpolation of the 
MC data. The quality of the DSM interpolation was usually superior to any single 
direct MC calculation apart from those with extremely high statistics. 
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Figure Captions 



Fig. 1 The minimal values of per degree of freedom for linear fits at fixed couplings 
A/g^ according to Eq. (11) for (|L|) (solid line) and Xv (short dashes). For 
the cumulant Qr (long dashes) Eq. (12) with uj = 1 was used. 

Fig. 2 The slopes f3/v (solid line) and ^/v (short dashes), and the constant cq = g"^ 
(long dashes) from the same fits as in Fig. 1. The dotted lines show the 
corresponding Ising model values, the dashed-dotted lines the error bars 

Fig. 3 The slopes of the linear fits to the logarithms of the derivatives of (|L|) (solid 
line), Xv (short dashes) and Qr (long dashes). The dotted lines show the 
corresponding Ising model values, the dashed-dotted lines the error bars 
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Figure 1 
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Figure 2 
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Figure 3 



